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TWO-DIMENSIONAL FINITE ELEMENT ANALYSIS 


OF RECTANGULAR PANEL WITH HOLE USING NICE /SPAR 

By 

1 2 2 
Zia Razzaq , Venkatesh Prasad , Siva Prasad Darbhamulla 

A. Panel Geometry and Loading 

A panel 30 in. long, 11,5 in, wide, with a 2,0 in, diameter hole at 
the center is analyzed. Since a two-dimensional analysis is conducted, 
the thickness of the panel is taken as unity. Figure 1 shows the panel 

and the applied compressive loading. Owing to the symmetry, it is 

sufficient to analyze only one-fourth of the panel with appropriate 
boundary conditions. 

B. Types of Discretizations 

Figures 2 through 8 show the various types of discretizations 

investigated using E41 quadrilateral elements of NICE/SPAR. Figures 9 

through 15 show the discretizations investigated using E31 triangular 
elements of NICE/SPAR. The triangular element discretization is obtained 
by simply adding a diagonal to each of the quadrilaterals. An enlarged 
view of the discretization near the hole is also given in each of these 
figures. Increasing number of elements are considered in order to 
conduct a systematic convergence study. 
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C. Isotropic Panel 

The following material properties are adopted for the isotropic 
panel : 

E - 10,000 ksi 
l> = 0.3 

which correspond to those of a typical aluminum alloy. Tables 1 and 2 

present a summary of the results for the quadrilateral, and triangular 

isotropic elements. In these tables, the number of elements used for 

one-fourth of the panel, the total CPU time for each computer run, the 

largest principal compressive stress, cr > and the corresponding 

v max 

element number are given. Figure 16 gives the graphical representation 

of the relationships between nr and the number of finite elements used. 

''max 

Clearly, the triangular elements provide a better estimate of the maximum 
stress (whose theoretical limiting value is 3.0). Figure 17 shows the 
relationships between CPU time and the number of finite elements. The 
difference between the CPU time with 120 quadrilateral elements and 240 
triangular elements is not dramatic. The information in Figures 16 and 
17 leads to the conclusion that the triangular elements are more suitable 
than the quadrilateral ones, for the two-dimensional stress analysis 
problem considered here. 

D. Ortho tropic Panel 

The following material properties are adopted for the orthotropic 
panel: 

E x = 10,000 ksi 
E y = 1,000 ksi 
\) =0.3 

Tables 3 and 4 present a summary of the results for the quadrilateral and 


2 


triangular ortho tropic elements. 


Figures 18, and 19 are the 


corresponding rv , and CPU time versus the number of finite elements 
max 

used, respectively. In this study, it is assumed that the NICE/SPAR 

program automatically accounts for the necessary transformations of 

material properties and stresses from the global to local coordinates and 

vice versa. Thus, the validity of the results for the orthotropic panel 

depends on the correctness or otherwise of this assumption. Figure 18 

shows a significant difference between the <T values for the 

'■'max 

quadrilateral and triangular discretizations. Figure 19 exhibits the 

same general character for the orthotropic panel as seen earlier in 

Figure 17 for the isotropic panel. Figure 20 shows how the value of cr 

max 

is affected due to a variation in the E^/E^ ratio when 160 triangular 

elements are used. For E /E = 1.0, the /r value from the plot agrees 

" x max 

with that obtained earlier using the isotropic triangular 

discretization. For other E /E values in the range from 0.1 to 1.0, a 

y x 

smooth curve follows. However, the validity of this curve also depends 
on whether or not NICE/SPAR is handling the necessary transformations 
properly. 

On-Going Study 

A three-dimensional stress analysis near the hole as well as the use 
of substructuring and subsequent parallelization of computations are 
included in the present research activity. 

Publication 

A paper titled "Concurrent Processing for Nonlinear Analysis of 
Hollow Rectangular Structural Sections," by Siva P. Darbhamulla, Zia 
Razzaq, and Olaf 0. Storaasli, has been accepted for publication in 
Engineering with Computers: An International Journal for Computer-aided 

Mechanical and Structural Engineering , 1987. 
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Table 1. Summary of results for rectangular panel with hole 

using E41 (NICE/SPAR) quadrilateral isotropic elements 


Number of 
elements 

CPU Time 
(secs) 

Element No. 
with CT* max 

cr 

max 

(ksi) 

12 

65.5 

4 

1.35 

22 

73.9 

7 

1.97 

36 

88.1 

7 

1.98 

48 

104.2 

11 

2.03 

68 

126.7 

11 

2.03 

80 

146.3 

22 

2.11 

120 

248.4 

6d 

2.13 


\ 
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Table 2. Summary of results for rectangular panel with hole 
using E31 (NICE/SPAR) triangular isotropic elements 


Number of 
elements 

CPU Time 
(secs) 

Element No 
with ^max 

O' 

max 

(ksi) 

24 

69.9 

8 

1.76 

44 

80.8 

11 

2.38 

72 

96.1 

11 

2.39 

96 

120.2 

15 

2.49 

136 

142.9 

15 

2.50 

160 

164.1 

29 

2.67 

240 

285.7 

71 

2.74 






Table 3. Summary of results for rectangular panel with hole 

using E41 (NICE/SPAR) quadrilateral ortho tropic elements 


Number of 
elements 

CPU Time 
(secs) 

Element No. 
with 0~max 

^max 

(ksi) 

12 

68.6 

5 

1.47 

22 

78.4 

7 

1.60 

36 

QQ 7 

yv . i 

7 

1 

1 . 62 

48 

105.5 

13 

1.64 

68 

127.1 

13 

1.61 

80 

146.7 

25 

1.66 

120 

243.7 

72 

1.95 


E * 10,000 ksi ; E - 

*■ y 


NOTE: 


1,000 ksi 
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Table 4. Summary of results for rectangular panel with hole 

using E31 (NICE/SPAR) triangular orthotropic elements 




























Figure 4. Discretization with 36 quadrilateral elements 
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Figure 15. Discretization with 240 triangular elements 
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